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Abstract: We propose renormalizable models of new physics that can explain various 
anomalies observed in decays of S-mesons to electron and muon pairs. The new physics 
states couple to linear combinations of Standard Model fermions, yielding a pattern of 
flavour violation that gives a consistent fit to the gamut of flavour data. Accidental sym¬ 
metries prevent contributions to baryon- and lepton-number-violating processes, as well as 
enforcing a loop suppression of new physics contributions to flavour violating processes. 
Data require that the new flavour-breaking couplings are largely aligned with the Yukawa 
couplings of the SM and so we also explore patterns of flavour symmetry breaking giving 
rise to this structure. 
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1 Introduction 

Over the last few decades, the Standard Model (SM) has emerged victorious as the correct 
description of physics up to and somewhat beyond the weak scale. Nowhere is this more 
true than in the flavour sector, where (with the addition of neutrino masses) the Yukawa 
interactions give rise to a pattern of fermion masses, mixings, and CP violation whose 
consistency with myriad precise experimental observations is striking. Crucial to this 
consistency is the particular (even peculiar) structure of the SM, which leads to many subtle 
and delicate cancellations in flavour-changing and CP-violating processes. As examples at 
tree level (see, e.g. [1] for more details), the representation content of the SM fermions 
(which is such that any two equivalent irreducible representations under the unbroken 
SU{3)c X U{l)em are also equivalent under the full SU{3)c x SU{2)l x U{1)y) ensures the 
absence of flavour-changing neutral currents (FCNCs) mediated by vector bosons, while the 
presence of just a single Higgs doublet field ensures the absence of FCNCs mediated by the 
Higgs boson. At loop level, the unitarity of the CKM matrix, together with the fact that 
it is empirically observed to be close to being diagonal, leads to remarkable suppressions 
of flavour-changing processes via the GIM mechanism. 
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Key to all of these suppressions, of course, is our insistence that the SM be renormal- 
izable; indeed, if we relax this requirement, we quickly run into conflict with data, unless 
the UV cut-off of the theory is rather high. 

In the last year or so, cracks have begun to appear in the SM edihce, in the form of a 
variety of anomalies associated with S-meson decays. These include effects seen in angular 
observables [2-10] of the decay B —)■ [11], the observable Rk [12], and a series of 

branching ratios with muons in the final state {e.g. [13, 14]). Whilst it is surely far too soon 
to claim that the end of the SM is nigh, it is perhaps worthwhile to explore how the SM 
might be modified in such a way that these anomalies in the data can be accommodated. 
Of course, this must be done in such a way as not to spoil the various delicate cancellations 
that occur in the SM, since this would lead to gross contradictions with data elsewhere. 
As such, it makes sense to look for theory which is renormalizable (or has a cut-off which 
is much larger than the scale of new physics, which amounts to the same thing). Then we 
can at least hope that, just like in the SM, dangerous flavour-changing processes can be 
kept under control. 

A number of models explaining the anomalies have already been proposed. These can 
be sub-divided into models which generate the anomalies via tree- vs. loop-level corrections. 
Models in the first category include leptoquarks [15-23] and Z's [24-39]. Not all of these 
models are renormalizable. Ref. [17], for example, describes a well-motivated model in 
which the electroweak hierarchy problem is solved by Higgs compositeness and SM fermions 
are partially composite. A light leptoquark can arise naturally as a Goldstone boson [40]. 
But if the compositeness scale is sufficiently high, the renormalizable limit is approximately 
recovered. There is one recently-proposed model in the second category, which generates 
contributions that explain the B anomalies and the anomalous magnetic moment of the 
muon, via a Z' with loop-induced couplings to muons [41]. 

In this work, we wish to propose a second type of renormalizable model in the category 
of those generating the anomalies at loop level. In order to obtain the flavour structure 
that seems to be required (in the quark sector at least) in an automatic way, we insist 
that the new flavour-violating couplings be linear in the SM fermions. We then survey the 
various possibilities for the BSM fields and find that basic phenomenological considerations 
(such as the insistence on an accidental symmetry stabilising the proton), together with 
a minimality criterion, lead to just two possible models. The models both feature two 
new scalar fields and a single fermion field, which couple to linear combinations of the SM 
fermions via Yukawa interactions. The two models are very similar in their phenomenology 
and so we discuss only one in detail. We find that, for suitable values of the parameters, a 
satisfactory ht to the anomalies can be obtained. 

The fit to the data requires that the new flavour-violating couplings be strongly hierar¬ 
chical (at least in the quark sector; in the lepton sector there is more room to manoeuvre) 
and moreover largely aligned with the flavour breaking already present in the Yukawa cou¬ 
plings of the SM. Thus, even more than in the SM, the low-energy theory seems to be 
crying out for a fundamental theory of flavour. Rather than attempt to find an explicit 
theory of flavour that does the job, we content ourselves with showing the plausibility of 
the existence of such a theory, by exhibiting patterns of flavour symmetry breaking that 
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can give rise to the required structure. We find that a variety of symmetries are possible. 
Interestingly, they give rise to a pattern of couplings similar to that obtained in theories 
featuring partial compositeness. 

The models have another feature of interest. By construction, they yield explicit UV 
completions that generate the non-renormalizable flavour structure that was identified in 
[42] as a viable one for explaining the anomalies. The fact that such UV completions exist 
is a desirable thing, because the alignment of the flavour-violating new physics couplings 
in [42] is not preserved by the SM RG flow. Thus, in the absence of an unexplained hne 
tuning in the couplings, the scale of new physics completing the effective lagrangian in 
[42] should be light, in order that the picture makes sense. Our models provide explicit 
completions. 

The structure of the paper is as follows. In Section 2 we describe the construction and 
content of the models. In Section 3 we outline phenomenological implications of one of the 
models and find allowed parameter space regions to ht the B anomalies. Section 4 contains 
an analysis of possible flavour breaking patterns, and our conclusions are given in Section 

5. 

2 The Models 

Recent experimental data in flavour physics, in particular measurements by the LHCb 
collaboration [11, 12, 14], suggest possible effects of New Physics (NP) in semileptonic 
decays of R-mesons. In particular, considering the 6 —)• s quark flavour transitions, the 
most significant departures from the Standard Model predictions are observed in: (i) the 
so-called angular observable of the B —)• K* fjP decays [11] (ii) a series of branching 
ratios of B-decays with muons in the final states (like the recently measured B^ —)• 

[14]) (in) the observable Rk [12] defined as the ratio of branching ratios of R —)• 
and R Ke~^e~ in the low (lepton pair invariant mass squared) region. It is perhaps 
premature to claim evidence for NP. Indeed the discrepancies in (i) and (ii) could be due 
to underestimates of hadronic uncertainties^, while the discrepancy in (Hi), despite the 
observable being theoretically very clean in the SM, could simply be a statistical fluctuation. 
However, if we allow for an interpretation of these experimental results in terms of NP, it is 
quite remarkable that model-independent approaches based on higher-dimension operators 
[5, 7, 8, 15, 24, 45, 49-54] give a simple and consistent fit to the data. 

In general, we expect that NP will couple chirally to the matter fields; assuming a 
coupling purely to either right- or left-handed currents, the fits find that the operator 
bLl°‘sLjlLloiHL is preferred. Motivated by this, we seek renormalizable models that couple 
the NP to quark and lepton doublets. Moreover, for reasons that will become clear, we 
demand that the SM fermions appear linearly in the NP couplings. In what follows, the 
SM fermions are written as Ql, Ur, Dr, Lr, Er while the beyond SM (BSM) fermions and 
scalars are denoted by Tj and <I>j. The Higgs doublet is denoted by H and transforms as 
(1, 2, 2 ). The possible linear interactions of the new BSM fields with the SM doublets {Q\ 

^For discussions on the size of the hadronic uncertainties on these observables we refer the reader to 
e.g. [6, 8, 10, 43-49]. 
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and -L^) can be classified in the following way: 


af + af + h.c.; 

(2.1) 

or + af+ h.c.; 

(2.2) 

or Q;fTgQ/d>-kaf^£L/4>-kh.c.; 

(2.3) 

or al + a\ + h.c. 

(2.4) 


Notice that one new state appears simultaneously in both the quark and lepton interactions 
(this is needed in order to draw a diagram contributing to 6 —)■ s^u/r). In the first two 
cases the common mediator is 'k while in the last two the common mediator is <I>. There 
are infinitely many combinations with suitable gauge quantum numbers to yield these 
interactions. In order to reduce the possibilities for the SM quantum numbers of the new 
states, we will impose conditions on them based on the following considerations: 

(a) Accidental protection 

The irreducible representations (‘irreps’) of new states should be such that other renor- 
malizable couplings with SM fermions beyond the Yukawa terms above are forbidden by 
gauge invariance. This guarantees that the definitions of Baryon and Lepton number 
can be extended to include the new sector, and that they will then remain accidental 
symmetries of the model. It also prevents other sources of flavour breaking. This crite¬ 
rion immediately implies that none of the new states can be in the same representations 
under the SM gauge and Lorentz groups as any SM field. 

(b) Scalar couplings with the Higgs doublet 

New scalars will always have quartic interactions with the Higgs doublet. In particular, 
terms of the form H) and are never forbidden by gauge 

symmetry (though the second one is absent if d> is a SU{2)l singlet). These terms 
are phenomenologically viable, but other quartic and trilinear interactions of a scalar 
with the Higgs, like ^HH, or ^HHH, could give rise to a violation of the custodial 
symmetry at the tree level and/or could modify the observed Higgs phenomenology. 
To be safe from these unwanted effects, we choose the quantum numbers of the new 
scalars such that these dangerous interactions are prohibited. 

(c) Direct searches, coloured particles 

New particles in the loops will need to be rather light to create a measurable effect in B 
decays, so it is convenient to choose quantum numbers such that their masses are less 
constrained by direct searches. The quark interaction requires at least one state that 
transforms non-trivially under the colour group. A coloured scalar will have weaker 
bounds on its mass than a coloured fermion, since in the latter case the production 
cross section is higher for a fixed gauge quantum number. Selecting a scalar to be the 
only new coloured particle leads us to consider only the first two cases in the above 
list of new interactions, namely the ones with a single fermion mediator 'k and two 
different scalars <k|j and $£. 
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Field 

SU{3)c X SU{2)l X U{1)y 

Interactions 


(1,1,0) 

or Wli 


(1,2,^) 

^iU^nQi 


(1,3,0) 

H^H 


(1,3,1) 

4>£ HtiLt 


(1,4,^) 

H^H^H 


(1,4, i) 


4- 

(1,1,0) 



(1,2,^) 



(1,3,0) 



(1,3,1) 



Table 1: Irreducible representations (with d < 5) of multiplets that contain a colourless 
and neutral particle but which are rejected because they could give rise to unwanted renor- 
malizable interactions, as listed in the last column. Irreps with negative hypercharges are 
related to these ones by charge conjugation. 

(d) Direct searches, BSM Lightest Particle (LP) 

The Yukawa interactions above are manifestly invariant under a U{\) transformation 
that acts non trivially only on the BSM states. This symmetry is respected by the 
gauge-kinetic terms of the new states too. We look for irreps such that this trans¬ 
formation is an accidental symmetry of the whole renormalizable model. This has 
the advantage that all NP flavour-violating processes are loop suppressed. But it also 
implies that the lightest NP state is stable; in order to evade strong constraints on 
coloured and/or electrically charged stable states coming from colliders and cosmol¬ 
ogy we look for a colourless irrep containing a neutral particle. The gauge quantum 
numbers are (l,re, ?/) under SU{3)c x SU{2)l x U{1)y and requiring an uncharged 
lightest state gives the restriction that —n/2 <y<n/2wyis integer if n is odd and 
is half-integer if n is even. 

(e) SU{2)l Minimality 

Finally we apply a minimality criterion on the dimensionality of the SU(2) irreps. This 
criterion is harder to justify on phenomenological grounds, though larger representa¬ 
tions will, of course, lead to a Landau pole at lower energies. For the sake of argument, 
we require that all new irreps have dimensionality fewer than 5. 

Let us now identify the irreps that could contain the LP. From (d) and (e) we obtain 
a hnite list of candidates. Most of them are excluded because of the presence of renormal¬ 
izable interactions that violate conditions (a) and/or (b). The excluded cases are listed 
in Tab. 1. We are left with four cases, each of which has a single fermion, tk, with SM 
quantum numbers (1,4, ±1/2) or (1,4, ±3/2). However, radiative corrections split the val¬ 
ues of the particle masses in the multiplet and it turns out that, for the quantum number 
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Field 

SU{3)c X SU{2)l X U{1)y 

Uil)B^ X U{1)l' X [/(l)x 

Ql 

(3,2,i) 

O 

O 

Ur 

(3,1,1) 

O 

O 

Dr 

(3,1,-1) 

(1,0,0) 

Ll 

(1,2,-1) 

(0,1,0) 

Er 

(1,1,-1) 

(0,1,0) 

4>h 

(1,2,1) 

(0,0,0) 

4' 

(1,4,-i) 

(0,0,1) 


(3,3,1) 

(-1,0,1) 

4>£ 

(1,3,2) 

(0,-l,l) 


Table 2: Quantum numbers of the Standard Model fields and new fields under the SM 
gauge symmetry (second column), and under the accidental global symmetries of the theory 
(third column). 

(1,4, ±1/2), the LP is not the neutral one. We conclude that, since we are demanding a 
neutral LP, the LP can only be contained in the fermion field 'L with quantum numbers 
(1,4, ±|). Imposing condition (e) on the field 4>q we are left with just two models: 

• Model A. 4' ~ (1,4, ±|), 4>g ~ (3,3,1), CD, r\j (1,3, 2) with Yukawa interactions as in 

( 2 . 1 ): 

aj^Qi^ + af^Li^e + h.c. (2.5) 

• Model B. 4^ ~ (1,4, —|), <I>q ~ (3,3,—|), ~ (1,3,2) with Yukawa interactions as 

in (2.2): 

af + af ± h.c. (2.6) 

The two models have very similar implications for the phenomenology that we are interested 
in here. Henceforth, we discuss only Model A. 

The quantum numbers of the SM and NP fields under the gauge and global symmetries 
(to be discussed below) are summarised in Tab. 2 and the most general renormalizable 
lagrangian is given by 


C — CsM ± ± Tvjr ± Tyuk, (2-7) 

u = ± - V{^h, ^g, ^e), ( 2 . 8 ) 

= i^L)^7^4' - M,k^4', (2.9) 

Ain = a^i^RQi<^q + ai^RLi<^i + aT~^L^R^l + atVL'i>R<^l. (2.10) 


See Appendix A for the explicit decompositions of the operators in terms of components 
of the SU(2)l multiplets. Let us now analyse the accidental global symmetries of this 
lagrangian. Before considering the breaking coming from C\\^ it is easy to show that the 
Lagrangian is invariant under a global 17(1)^. Indeed, the SM alone has accidental global 
symmetry U{1)b x t7(l)e x U{1)^ x U{l)r, while the gauge kinetic terms of the new BSM 


- 6 - 








fields have global symmetry C/(l)^ x [7(1)$^ x 17(1)$^. Moreover, it is easy to prove that 
the most general renormalizable scalar potential is invariant under U{iy. 

Now consider the effect of C\{^. For a generic choice of the couplings and there 
is always an unbroken 17(1)^ = U{1)b’ x 17(1)l' x 17(1)x) defined as follows. Under the 
U{1)b' the SM fields have their usual baryon number while has charge -1/3. Similarly, 
under the 17(1)/,/ the SM fields have their usual lepton number while has charge -1. 
Finally the SM fields are uncharged under 17(1);,^, while the BSM fields have charge unity^. 

Thus, the model retains analogues of the accidental baryon and lepton number sym¬ 
metries of the SM, which suffice to stabilize the proton and to prevent contributions to nu¬ 
merous unobserved lepton- and baryon-number violating processes. Moreover, the model 
features an additional accidental 17(1)^ symmetry, under which SM fields are uncharged. 
An immediate consequence of this is that all NP-generated processes involving only SM 
particles in the initial and final states are loop-suppressed. This is certainly an advantage 
from the point of view of the vast majority of flavour-violating observables, where no de¬ 
viation from the SM is observed. It might be regarded as a disadvantage from the point of 
the view of the B-physics anomalies, where a sizable NP effect is needed. But this is offset 
somewhat by the desirable structure of linear NP flavour violation that results. We shall 
see in the sequel that the anomalies can be reproduced even for values of the NP couplings 
that are of order unity or smaller. 

3 Phenomenological Analysis 

In this section we discuss the phenomenology of Model A. In an obvious notation, we 
denote the masses of the new states as Mij/, Mg, and M^. In a basis where the left-handed 
quark doublet is defined as = {VcKMu\,d^j^)'^, the minimal set of couplings a that are 

needed to fit the b —)• sii anomalies are Og, and (ie. couplings involving b, s and 
//). To begin with, we will assume that only these couplings are non-zero and investigate 
the processes induced. In this section we collect relevant formula on indirect searches, and 
investigate direct production bounds. In subsection 3.3 we use this information to find 
allowed parameter space regions. In section 4 we will discuss relaxing the assumptions on 
the couplings and propose more motivated flavour structures. 

3.1 Indirect searches 

As described above, the accidental global 17(l)x symmetry under which the new particles 
are charged implies that contributions to processes containing only SM particles in the 
initial and final states are only induced at loop level. Here we investigate the size of these 
contributions. 

The relevant processes, given the assumption on couplings described above, are b —)• 

processes, Bg mixing, b —)• sy, and the anomalous magnetic moment of the muon. 

^This symmetry makes our neutral and colourless LP stable. Hence, the LP is a potential dark matter 
(DM) candidate. However, even if its mass and couplings could be fixed in order to reproduce the right 
relic abundance, data from direct detection experiments exclude it as DM, because of its large coupling to 
the Z-boson. A potentially large relic density could be problematic, though non-renormalizable operators 
are expected to break the accidental U{l)x at very high scales, triggering its decay into SM particles. 
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Figure 1: Diagram contributing to 6 —)• sfj,^ 


3.1.1 Semileptonic four-fermion operators 

The process b —)> sii, important for the LHCb B meson anomalies, is induced at loop 
level by the diagram in Fig. 1^. The SU{2)l structure of the NP-induced semileptonic 
four-fermion interaction can be derived from the discussion in Appendix A, using the 
lagrangian (Eqn. A. 6) written explicitly in terms of SU{2)l components. The resulting 
effective NP lagrangian is 


(3.1) 

with Xq = jA and = jA. The loop function K{xq, x^) can be obtained by the following 
definitions; 




K{xq,X() 


Ml 


64 TT^ 


{q\ 


K{x) 

K{x,y) 


1 — X + x"^ log X 
{x — 1)2 
K{x) - K{y) 
x-y 


The effective hamiltonian relevant to 6 —?■ sU transitions is 


^ Cl{y) Oliyi) , (3.2) 

where are a basis of SU(3)c x U(l)Q-invariant dimension-six operators giving rise to the 
flavour-changing transition. The superscript ^ denotes the lepton flavour in the final state 
{I G {e,/i, r}), and the important operators for our process, Of, are given in a standard 
basis by 

^ ^ (sjaPL(R)b) , (3.3) 

^ {siaPLiR)^) 

^There are also Z and photon penguin diagrams which contribute, with a NP loop connecting the quarks 
and joining to the leptons via a ^/7 propagator. These penguin diagrams are discussed in Appendix B and 
are found to be very suppressed relative to both the SM contribution and the diagram in Fig. 1, and hence 
are neglected here. 
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Figure 2: Diagrams contributing to Bg mixing 


Comparing equations 3.1 and 3.2 we find the NP contribution to the Wilson coefficients 
relevant to 6 —)■ is 


'-^9 ~ '-"10 


ot \ ^ 7 K{xq,xe) 

~f7r^tsnbj- r-ryc-'?. 172 


VV2 


dvr J 57671^ 


Ml 


Q* q 
O^o O^Q 


ai 


(3.4) 


The most recent best fit ranges on this combination of Wilson coefficients are taken from 
[49] and are given by 


C^NP ^ g [-0.71, -0.35] (at la), 

C^NP ^ _ci^^P g [-0.91, -0.18] (at 2u). 


(3.5) 

(3.6) 


3.1.2 Four-quark operators 

Interactions between four quarks are induced at loop level by diagrams like those in Fig. 
2. These interactions can lead to meson mixing; in particular, if the process b —)• sfifi 
is present, then inevitably Bg mixing must also be induced. This process can therefore 
introduce important constraints on the masses and couplings of the new particles. The 
four quark effective operator induced by the NP is 


^eff D 


K'{xq) af a^^aVtO-n 
Ml 128 7r2 


{QTi,Q 1) + ^ • {QlImQI) 


(3.7) 

where K'{x) is the hrst derivative of K{x). The SU{2)l structure of the effective operator 
is similar to that of Eqn. 3.1 and can again be derived from the discussion in Appendix A. 
Projecting the quark doublet along the down components we find that for Bg mixing the 
relevant operator is 


Ceff A m2^ («r«3)^ {sLl''bL){sL-ft,bL) + h.c.. (3.8) 


The Wilson coefficient is easily extracted at high energy /r = A where the BSM particles 
are dynamical fields. We fix A = 1 TeV in what follows. At this energy we have 


C?"(A) 


7 K'jXq) 
5767r2 Ml 




Uo a 


3J 


(3.9) 
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Figure 3: Diagram contributing to 6 —)■ sj. The photon is attached in all possible ways. 


In order to place bounds on the parameters of our model, we take into account QCD 
effects using the results and procedure of [55]. Using the anomalous dimension of this work 
we found that the running of Wilson coefficient from the scale of the New Physics (A) to 
the scale of the process (rufe) is given by = rjvLLCi^i^) with tjvll = 0.78. For 

the evaluation of the relevant matrix element we used the lattice result of [56] . These lead 
to a constraint (at 95% conhdence level) on the coefficient 


C^(A) < 


1.8 X 10-5 

TeV^ 


which translates into 


7 K'jXg) 
5767r2 M| 


al*alf < 1.8 X 10-5 


1 


1 TeV 


2 


(3.10) 


(3.11) 


Thus the measurement of Bg mixing produces a bound on the hadronic couplings 
involved in the b s^fj, process, viz. and Ug. The model can hence accommodate both 
this bound and the b —)> s/i/U data if the muonic coupling a\ is sufficiently large. In this 
respect the model is similar to Z' models — the couplings involved factorize into leptonic 
couplings and hadronic couplings which can be set independently. This factorization does 
not occur in leptoquark models. 


3.1.3 6 —)• S7 

The radiative process 5 —>• sy will also be induced by the diagram in Fig. 3. The couplings 
involved are the same as those for Bg mixing. However, the amplitudes will scale differently 
with the parameters Uq and Mq between the two processes. Constraints from b ^ s'j could 
therefore provide complementary information. 

At the mass of the b quark, the process 6 —?• sy is described by the following effective 
hamiltonian: 

-Hes = {V^gVth) [CrimhWmt) + C^(mb)0'(m^)] , 

where {saa 0 PR{L)b) . 
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Figure 4: Diagram contributing to the anomalous magnetic moment of the muon. The 
photon is attached in all possible ways. 


At the matching scale M, we get an additional contribution from the NP to the coef¬ 
ficient of the dipole operator; 


Gf, 


where Fi{x) is defined as 


— 1 Q* Q 
Q-o 0^3 

I 




(3.12) 


Fi{x) = 


[x^ — -|- 3x -|- 2 -|- 6x log x) 


(3.13) 


12(x- 1)4 

The 2fT allowed range for this parameter has been fitted recently in [49], giving 

Cf^^(mfc) G [-0.10,0.02] (at 2u). (3.14) 


3.1.4 Anomalous magnetic moment of the muon 

Although it is somewhat peripheral to our discussion, let us remark that loops of T and 4>£, 
as shown in Fig. 4, generate a 1-loop contribution to the magnetic moment of the muon, 
which may be able to resolve the long-standing experimental discrepancy therein [57]. The 
NP contribution is given by 


\al\‘^ Mf, f 2 1 \ 

=^4r4.,)+-u(.,-‘)). (3.15) 

which should be compared to the observed discrepancy [58] 

Ao/, = = (287 ± 80) x 10“^^ (3.16) 

As we will show in Section (3.3), it is possible to fit the anomalous magnetic moment in 
this model. However, it requires a large value of which is problematic, since it can lead 
to large corrections to electroweak precision observables at the Z-pole. 


3.1.5 h —)■ svv processes 

Contributions to B ^ KvV and B —?• K*fV are expected in the model, due to a diagram 
similar to Fig. 1 with the muons replaced with muon neutrinos (as well as Z penguin 
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diagrams — see the comment in section 3.1.1). A detailed analysis of NP contributions 
to this process is given in [59], and we use their results here. Current measurements give 
bounds on the ratio of total (NP+SM) to SM branching ratios to be 


_ BR{B KW) 

^ ^ BR{B Kuu)sm 


(3.17) 


TZk* 


BR{B K*Vv) 
BR{B —^ K*vv)sm 


< 4.4 


(3.18) 


at 90% confidence level. An expression for 77^(*) in the presence of NP with couplings to 
left-handed SM fermions is 


TIk = TIk* 



\ci? 


\CSM\2 


(3.19) 


where is the coefficient of the operator 

Ol = (si^Pib) (Pn'^(l - 75)^i) (3.20) 

in the effective Hamiltonian. The SM Wilson coefficient is known to be quite accu¬ 

rately = —7.65. In the case of our model, 


C 


NP,fi _ f 


/ vr \ 

/ 4K(xg,xe) 


yV2GFaVaV^*J 

V 9 Ml 



(3.21) 


where the definitions of K{x,y), Xq and xi are as in Sec. 3.1.1. Thus bounds exist on the 
parameters of the model due to 6 —>• suv processes: 


- 24.9 < < 30.0, (3.22) 

-55A<K{xq,xe)alaf\ai\^ < 66.8, (3.23) 

at 90% conhdence level. However, if we assume that couplings to muons are the dominant 
leptonic coupling in the model, then we find the relation 

Cr = (3.24) 

Therefore, for values of the Wilson coefficients required to fit the b —)• s/r/r anomalies (Eqn. 
3.5), the NP contributions to the branching ratios for the b —)• suv processes are well below 
the bounds, adding approximately 5% to the SM values. 


3.2 Direct Searches 

The particles of the three new multiplets, 4>g, <I>£, and T, will be directly produced at the 
LHC if their masses are within kinematic reach. In this subsection, we outline current 
limits on their masses from direct searches, and identify promising channels to search 
for them. It will be convenient to label the SU{2)l components of the multiplets by 
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a superscript denoting their respective electric charges; the full list of NP particles is, 
therefore, ^'=(^'+3,^+2,^+i,^r0). 

The collider phenomenology of the new particles will depend on the mass spectrum. 
As before, we assume that < Mi^Mq, since we require the LP to be neutral. We 

further assume that Mi < Mq, since this minimises contributions to Bg mixing (it also 
maximises contributions to the muonic g — 2). The three multiplets can, generically, be 
well-separated in mass, but within each SU{2)l multiplet there may also be significant 
mass splittings. For the scalar multiplets, there are tree level mass splittings due to the 
presence of direct couplings with the Higgs; for the fermion multiplet, there are only small 
radiative mass splittings between the components. In the limit that the common mass is 
much larger than the electroweak scale v, the radiative mass splitting between the different 
charge eigenstates is [60, 61] 

^"irad = mQ+i - niQ ^ 166MeV ( 1 -h 2Q H-) , (3.25) 

V cos Oyi) 

a formula which holds for both scalars and fermions. According to eq. 3.25 the lightest 
particle within the fermion multiplet will be uncharged, as desired. 

As the 'k fermion multiplet has the lowest common mass of all the new states, and due 
to the U{l)y^ symmetry within the new sector, the lightest state within the multiplet will 
be stable. The small radiative mass splitting means that heavier fermion components will 
decay to the lightest (neutral) component by emission of one or more soft charged pions 
or leptons, which will not be energetic enough to be reconstructed in the detector. Thus 
if any T particle is produced at the LHC, it will appear as missing transverse momentum, 
similarly to the Wino-like dark matter described in [62]'^. We thus neglect henceforth 
the soft undetectable pions or leptons emitted in the decays of heavier components of 'k. 
Therefore, and particles (being heavier than the T states) will effectively decay to a 
SM particle plus missing transverse energy. Furthermore, due to the U{\)^ symmetry, NP 
particles will always be pair-produced at the LHC. This means that searches for (R-parity 
conserving) supersymmetry should be sensitive to d>g and $£. We will now discuss each of 
the fields <kg, and T in turn. 

The fermions can be pair produced via a photon or a WjZ, or through the decay 
of ^q^i . By the arguments above, they will always behave as uncharged weakly interacting 
particles. Limits can be set on these from mono-x searches, from constraints on the invisible 
width of the Z boson, and from LEP searches for charginos that are almost degenerate 
in mass with the neutralinos. A detailed analysis of all of these has been performed in 
[63]; the last-mentioned has been found to be the most constraining, implying a bound of 

> 90 GeV. 

Each component (<k^^^^, <k^^^^) of the coloured scalar multiplet will be 

strongly pair produced and will decay with a similar signature to that of a squark; i.e. 

^Note that our setup is subtly different from that described in [62]. There, strong constraints can be 
put on the mass of the new fermion multiplet from disappearing tracks searches, since the lifetime of a 
charged fermion decaying to the neutral fermion can be long enough to create a disappearing track in the 
detector. Here, these searches are not constraining because the lifetime is too short for a track to be visible. 
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Figure 5: Diagram to illustrate the possibility of pair produced particles decaying to 
two different quarks. Since the pion produced when the decays to is too soft to be 
detected, the decay of the appears similar to that of a (anti-)sbottom, whereas the 
1 /3 

decay of the appears similar to that of a stop. 

to a quark and (either directly or via emission of soft pions or leptons) the stable neutral 
component of the fermion Note that there is the possibility that two d>q particles pro¬ 
duced each decay to a different flavour of quark. An example decay chain displaying this 
property is illustrated in Fig. 5. This complicates the re-interpretation of SUSY search 
limits. 

The flavour considerations discussed earlier only constrain the product of couplings 
Ogolj without constraining their quotient. The strongest constraints from existing LHC 
searches will hold for situations where one of the couplings or is much larger than 
the other, so that the branching ratio to a particular generation dominates. If Og ;§> 
the <I>g particles will decay like sbottoms and stops, whereas if Og <C they will decay 
like second-generation squarks. The branching ratio to an up-type quark as opposed to a 
down-type within a particular generation is determined by the SU (2 )l structure. 

We focus on the case Og 3> since this is motivated by flavour considerations, as 
explained in the next section. In this limit, one can show that the total branching ratio 
times cross-section for a pair of <I>q particles to be produced and to decay to a pair of tops 
is 7/8 of that for direct stop pair production. The most recent limits on direct stop pair 
production are given in [64, 65]. As a conservative estimate, given that we have the limit 

> 90 GeV, we can take the limits on direcly pair-produced stops decaying to tops and 
90 GeV neutralinos to apply to our 4>q. This gives a limit of Mq > 750 GeV. Likewise, the 
total branching ratio times cross-section for a pair of particles to be produced and to 
decay to a pair of b quarks is 7/8 of that for direct sbottom pair production. Latest limits 
on direct sbottom pair production are given in [66], and again taking these limits to apply 
to our particles, we find that Mq > 720 GeV for a T of mass 90 GeV. 
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Figure 6: Exclusion plot for direct pair production of sleptons decaying to a lepton and 
a neutralino, taken from [67]. The region under the blue line shows the exclusion on our 
model found by reinterpreting the exclusion plot in terms of direct pair production of $^8 
decaying to a lepton and a T. For our model the x-axis should be taken to mean the mass 
of the Mg, while the y-axis means the mass of the T, Mij,. The dotted part of the blue 
line is extrapolated. 


The SU{2)l components of the scalar ^g will, if they have only muonic 

couplings, always decay to either a muon or a muon neutrino, together with a T particle. 
So they will sometimes decay in the same way as a smuon in supersymmetry. However, the 
production cross-sections and branching ratios will differ. Results of recent LHC slepton 
searches are given in [67, 68 ]. These rule out left-handed smuons, pair produced directly via 
a VF/Z/y and decaying to a muon and a neutralino, up to a maximum mass of roughly 300 
GeV (for a massless neutralino). We used Feynrules [69, 70] and Madgraph5_aMC@NLO 
[71] to calculate electroweak (EW) pair production cross-sections of the ^g particles, and 
then multiplied by the branching ratios in order to reinterpret the limits on cross-section 
given in [67]. The CMS limit plot, with our model superimposed in blue, is shown in 
figure 6 . If the mass of the T particle is greater than about 150 GeV, there are no bounds 
on the mass of the ^g (other than the assumption that its mass is greater than that of T). 

3.3 Parameter space plots 

In this subsection we show allowed regions in the parameter space of the model considering 
the observables described above; b —)• sii, B —)• V^y, Bg meson mixing and the anomalous 
magnetic moment of the muon. The relevant parameters entering the expressions of these 
observables are the masses of the new states and Mg) as well as the coupling to 

muons 02 the combination 02 * 0 : 3 . Without loss of generality, we can re-define the 
BSM fields to make these parameters real. 

In order to fit the b —)> sii anomalies without being in disagreement with the measured 
Bg mixing rate, the muonic coupling 02 must be rather large. In order to have an idea of 
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the typical values of the parameters needed in our model, in Fig. 7 we show parameter space 
regions assuming that = 1-2 while parametrizing the masses in terms of one single scale 
M assuming the following hierarchy M-q, = M, = M + 200 GeV, Mq = M + 700 GeV. 
In this way we are left with two parameters only (M and The B —?■ Xg'y allowed 

region is not shown because it yields weaker constraints than Bg mixing does. For this 
hierarchy of masses, the only relevant direct production constraint is the bound on the 
mass of 'k, Mij, > 90 GeV. There is an overlap between the allowed Bg mixing region and 
the Icj preferred region for the b —)■ s£i measurements — so with these parameters, the 
model can fit the b —)• sii anomalies. The value of can be further lowered to be < 1, 
in this case the values of , Mi and Mq are close to present bounds coming from direct 
searches. For example we verified that a fit to the data with ~ 0.8 could be achieved 
when Mil, = 150 GeV, Mi = 200 GeV and Mq = 800 GeV. 

However, if we also wish to fit the anomalous magnetic moment of the muon, the muonic 
coupling a 2 must be larger. We show in Fig. 8 the relevant parameter space regions when 
this coupling is set to = 2.5, with the same hierarchy of masses as before. If we want to 
take this explanation of the {g — 2 )^ anomaly seriously, then we should consider possible 
bounds from the shift of the EW gauge couplings and W ~^(see also the 

discussion in section IV of [41]). The corrections are non-universal and so a global fit to 
EW data is required to establish the precise constraints on the couplings. Though such a 
fit is beyond the scope of our work, naive arguments suggest that 0 ( 1 ) values of are not 
problematic. 

4 Flavour symmetries 

In this section we establish a possible connection between the flavour violation present in 
the SM and in the NP sector. 

In the SM and in the limit of vanishing Yukawa couplings, the largest group of unitary 
field transformations that commutes with the gauge group and leaves invariant the kinetic 
terms is U{3)^ x U{1)h- Adopting notation similar to [72] we can decompose this group 
in the following way; 

Gk = SU{3fq X SU{3)j X U{1)b x U{1)l x 17(l)y x U{1)pq x U{l)En x U{1)h, 


with 


SU{3fq = SU{3)q, X SU{3)u^ x SU{3)d^ (4.1) 

517(3)1 = SUi3)L, X SU{3)e^. (4.2) 

The 17(1) factors can be identified with the baryon (B) and lepton (L) numbers, the 
hypercharge (Y), a transformation (PQ) acting non trivially and in the same way only on 
Dp and Ep, and finally a universal rotation for the fields Ep and a 77(1) global symmetry 
associated to the Higgs doublet. 

We would like now to make connections with the flavour structure of the SM and 
the possible effects coming from NP. In order to do that a first step is to identify (i) 
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Figure 7: Parameter space plot for = 1-2, and with the masses of the three fields given 
by Mq, = M, = M + 200 GeV, Mg = M + 700 GeV. For this value of a^, it is not 
possible to explain the anomalous magnetic moment of the muon whilst fitting the other 
constraints. 


a flavour symmetry and (ii) a set of irreducible symmetry-breaking terms. The flavour 
symmetry group Qp V Qk has to be broken in order to reproduce the observed pattern of 
fermion masses and mixing. In order to do that a set of symmetry-breaking spurions are 
introduced to formally restore the symmetry Qp- 

We will now consider 3 explicit examples and we will focus on the quark sector. 

1. gp = 17(3)3 

This is the case of Minimal Flavour Violation [72]. The spurion fields are the three 
Yukawa couplings 

Vf/~ (3,3,1) Yd ~ (3,1,3), (4.3) 

where the quantum numbers are specified with respect to the direct product of groups 
SU{?,)q, X SU{2,)u^ X SU{2,)d^. 

2. Qp = U{2fg 

This is the flavour symmetry of the quark sector if only the Yukawa couplings yt and 
yi) are non-vanishing. So to a good level this is an approximate symmetry of the SM. 
Recent works [73-75] considered the following set of irreducible spurions ; 

A„~(2,2,l), Ad ~ (2,1,2), V~(2,l,l), (4.4) 
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Figure 8: Parameter space plot for = 2.5, and with the masses of the three fields given 
by Mij = M, Ml = M + 200 GeV, Mq = M + 700 GeV. With this large value of there 
is an overlap between the regions that fit the B anomalies (in blue), and the anomalous 
magnetic moment of the muon (in green). 

where the quantum numbers are specified with respect to the direct product of groups 
SU{2)q^ X SU{2)u^ X SU{2)d^. 

3. Gf = U{lf 

This case mimics partial compositeness. The irreducible spurious are connected to 
the Yukawa couplings in the following way; 

{Yuh ~ efel {Ynh ~ ( 4 - 5 ) 

With these specific cases in mind we are now ready to discuss flavour violation induced 
by operators of the form a? and These operators break the 

flavour symmetry and in order to restore it we could assume that the vectors are again 
spurious with definite transformation rules under the flavour symmetry. We could now 
assume minimality of flavour violation in the following sense: the af can be expressed 
using the irreducible spurious used to construct the SM Yukawa couplings. Following this 
procedure we obtain the following results. 

1 . Gf = UiSfq 

To recover flavour invariance the have to transform in the following way; 

a'?~(3,l,l), (1,3,1), a‘'~(l,l,3). (4.6) 
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However, it can be proved using triality properties of the SU{3) irreps that tensor 
products of Yu,Y^,Y£),Y^ can never give rise to any of the . Furthermore, the 
structure of the Yukawas cannot be reproduced by combinations of the as, since this 
can only lead to Yukawa matrices of rank 1, meaning two of the up-type quarks and 
two of the down-type quarks would have zero mass. We therefore conclude that this 
structure cannot work. 

2 . Gf = 17 ( 2)3 

To recover flavour invariance the have to transform in the following way; 


(af,al)~(2,l,l), al~ (1,1,1), (4.7) 

(a)‘,af)~ (1,2,1), a^~ (1,1,1), (4.8) 

(af,a^)~ (1,1,2), af~ (1,1,1). (4.9) 

Atleading order in the number of spurion fields we have that 

= (4.10) 

(a)‘,a^), = a“(VtAJi, (4.11) 

= (4.12) 

while af’ = with , order one numbers and i = 1,2. Doing a spurion analysis 
with this setup, we find that in the basis in which = (ycKM'^Y^'^ for 1 = 1, 2 


(generation index) and Q\ = [e^^HL,bL)^, the coupling vector in the quark sector 
Oq is found to be 


ag = {yi 2 Vtd, yi 2 Vts, ys ), (4.13) 

where yi 2 and y^ are each (generally complex) numbers of 0(1). 

This is very similar to the case of Partial Compositeness, except that there is slightly 
less freedom in the couplings; the ratio of to af is exactly fixed (rather than fixed 
up to a factor of 0(1)). The natural size of agO^ is 0(|V)s|) ~ This ties in well 
with the allowed regions in Figs. 7 and 8. 

3. Gf = U{lf 

It is easy to show that link between the af and ef is simply given by af = cf ef 
where cf are order one numbers. 

Thus, a number of patterns of flavour symmetry breaking can give rise to the hierar¬ 
chical structure and alignment with the SM Yukawa couplings that is needed to explain 
the anomalies. In all cases, the pattern of couplings is similar to that arising in models of 
partial compositeness (see for example [76]). Although our phenomenological analysis in 
the previous section was done with more restrictions on the couplings, we can be sure that 
a partial compositeness-like flavour structure in the quark sector is phenomenologically 
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viable, since a full analysis of existing bounds on this setup has recently been performed 
in [17]. 

An analogous discussion could be repeated in the leptonic sector, starting from the 
global symmetry U{3)l^ x U{3)eh of the gauge kinetic terms associated with left- and 
right-handed leptons. However, a complete understanding of the leptonic flavour requires 
the knowledge of the mechanism that generates neutrino masses. Here we remain agnostic 
as to the possible flavour orientation of the spurions and a®. For recent works that 
consider possible links between the flavour violation in the neutrino sector and the physics 
of LFV transitions in 6 —)• sii', we refer the reader to [19, 42, 77-79]. 

5 Conclusions 

We have presented renormalizable extensions of the SM that can explain several anomalies 
observed in H-meson decays. Renormalizability (which amounts to the assumption that 
further NP is heavy and decoupled) allows us to introduce deviations from the SM, coming 
from NP, in a controlled way. This is almost a sine qua non, given that we observe just a 
handful of anomalies in data, while many thousands of other observations agree with the 
SM. 

We have surveyed the possible NP fields that allow for a coupling to linear combinations 
of left-handed SM fermions, since this generates, at one-loop, an operator of the form 

Coupled with the plausible assumption that the linear flavour- 
violating spurions, are roughly aligned with the Yukawa couplings of the SM, we end 
up with a good fit to the anomalies, without contradicting other data. 

As a spectacular example of the control that renormalizability brings, the models 
that we identify feature 3 accidental global symmetries, corresponding to conservation of 
(generalized) baryon and lepton numbers and to a ‘NP number’. The consequences of 
these accidental symmetries are manifold. Not only is the proton stabilized, but also all 
other baryon- and lepton-number violating processes {e.g. neutron-antineutron oscillations, 
IX —>• evv), many of which are strongly constrained, are forbidden automatically. The NP 
number leads to a generic suppression of NP flavour-violating processes, since these can 
only occur at loop level. Yet another advantage of the models is that NP can only couple 
to left-handed SM fermions at the renormalizable level, meaning that contributions to 
processes requiring a helicity flip, such as /r —)• ey, are further suppressed. 

The accidental symmetries of the models, while sufficient to prevent many dangerous 
processes, are quite different from the accidental symmetries of the SM, namely baryon 
and individual lepton family numbers. This is a crucial feature, since it allows us to have 
large violations of lepton universality. This is precisely what is needed to fit the anomalies. 

The models are not panaceas, in that there is a further anomaly in B physics that 
cannot be explained, arising in decays to [80-82]. But it seems hard to explain this 

anomaly in any NP model, for the simple reason that the SM contribution, with which 
it needs to be comparable, is so large (being a tree-level effect with minimal Cabibbo 
suppression)^. 

^However it is not impossible, see for example [21, 83, 84]. 
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The main weakness of the models is arguably that they require a rather large value of 
the coupling in order to explain the anomalies. While this coupling can be rather smaller 
than other couplings in the flavour sector [i.e. the top quark Yukawa coupling), some 
readers may be alarmed that such a large coupling should appear in the light lepton sector, 
where (at least in the SM) all other couplings are small. It is important to note, however, 
that not only does this coupling not cause phenomenological problems per se,® but also that, 
provided that it is suitably aligned, it does not lead to large flavour violations in the light 
leptons via renormalization group flow. This follows immediately from the observation that 
there exists a basis in which the SM leptonic Yukawa couplings are diagonal. Nevertheless, 
the necessary alignment is aesthetically disturbing; we have shown that it is plausible from 
the point of view of flavour symmetries, but it would be nice to have an explicit model of 
flavour in which it is realised dynamically. 


Acknowledgments 

This work has been partially supported by the Galileo Galilei Institute for Theoretical 
Physics, STFC grant ST/L000385/1, and King’s Gollege, Gambridge. We thank members 
of the Gambridge SUSY Working Group for discussions. 


A SU(2)i, decompositions 


By denoting the generators in the fundamental representation of SU(2)i as T“ = cr“/2 
(with ct“ being the Pauli matrices and a = 1, 2, 3), we define their action on the (2j + 1)- 
dimensional completely symmetric tensor (UG 2 , ■ ■ ■ ,i 2 j = 1, 2) as 


<5“(Xji*2...*2j ) — '^hk Xki2...i2j + T^^k Xiik...i2j + • • • + Xiiii...k 


(AT) 


In general, we arrive at the following embedding of the properly normalized electric charge 
eigenstates: 


X11...1 = 


X11...2 = 


X12...2 = 


X22...2 = 


\/ ■B2i,0 
1 

\/B2j,l 


Q 

Q-1 


(A.2) 


\/ ^2j, 


2j-l 


=x 


Q-2i+i 


Q-2j 


where the superscripts denote the electric charge of the field, Bn,k is the binomial factor 

D _ n\ 

On,k — k\(n-k)\ 


and the normalization of the states is such that 


X"l*2-*2^Xili2...i2i = IX^P + IX^ T + --- + |X + ^1^' 


(A.3) 


In the following, we provide the SU(2)i decomposition for the BSM fields ('I’,<kg and 
<1>£ ) introduced in the Model A (2.5): 

®In fact, as shown in [85], one can even put 0(1) couplings among the light quark generations without 
necessarily getting into trouble. 
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'I/ =(1,4,§) 

4'iii = 4'^ 

4^112 = 

4/122 = ^4/1 
4-222 = 4-° 


% =(3,3,I) 

(4),)n = 

(4>,)i2 = 


(4>g)22 — 4>q 


1/3 


4>£ =(1,3,2) 

(4>^)ii = 4-3 
(4 >£)i2 = ^4>| 
(4>£)22 = 4>^ 


The relevant linear interactions (2.10) introduced in Model A can be rewritten in the 
following way 


/:iin = ^rQI% + al + h.c. 


(A.4) 


= + h.c., (A.5) 

where A;i,/c2,^3 = {1,2} are SU{2)r fundamental indices. More explicitly we get 


Ain = al ( 4-^^4>^/3 ^ 2 ^^ 2 ^ 4/3 j ^ 


1:^-1: 


+ a. 


yi + \ll ^r^T + dl + h.c. 


I ,Ti 3^3 I , /2 ,J, 2^2 I , / ^ ,T, 1/F>1 


+ “i 1 4-^ 4>£ + y - 4/^ 4)^ + Y - 4-^ 1 1 / 1 + h.c. 


+ 


1 .T. 2^3 I . /2 ,y, 1^2 I , T ,0 ^1 


- 4/^ 4>£ + Y - 4/^ 4>£ + 4/^$^ 1 el + h.c. 


(A.6) 

(A.7) 

(A.8) 

(A.9) 


B Photon and Z-boson mediated contributions 


An explicit calculation of the photon penguin diagrams leads to the following (lepton 
flavour universal) contribution to the Wilson coefficient of the leptonic vector current at 
low energy; 


^9" = 


^^V*V — 
V2 *'471 


-1 


aT a% a 2 


Ml 


47r27A 


f'yix 


-l\ 


(B.l) 


where 


f'yix) = 


2(38- 79x + 47x^) 
9(1 — x)3 


4(4 — 6x + 3x3) 
3(1 -x)4 


logx, 


(B.2) 


and where, as before, Xg = This function is normalised in such a way that /^-(l) = 1. 

In our model, the contribution of the Z-boson mediated penguin diagram is suppressed 
compared to that of the photon mediated one. Indeed the Z-boson exchange is enhanced 
only in diagrams containing a source of explicit SU{2)r breaking^; such diagrams are not 
present in our case. When there is no explicit SU (2 )r breaking, the contribution from the 


^This general argument has been given in the context of Z-penguin contributions in the MSSM in [86] 
and in [87]. 
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Z-boson penguin diagram is suppressed by a factor ^ ~ 3 x 10“^ compared to the photon 
one; this factor is given simply by the ratio of the propagators in the two cases. 

Neglecting the Z-boson contribution, we now quantitatively show that the photon 
contribution is suppressed at the percent level when compared to the one in Eq. (3.4). 
Taking the ratio of the Wilson coefficients generated through a photon penguin and the 
NP box diagram respectively, we get 








Q 


PlNP 


WfhK{xq,Xi) \ai\ 


= 3.5 X 10“^- 


Xq f'^/jXq^) 

K{xq,xi) lagi 


(B.3) 


Taking as a reference the benchmark defined in Figure 8, namely = 2.5, Mq/ = M, = 
M + 200 GeV, Mq = M + 700 GeV we find that 


'-'9 


a 




< 3.9 X 10 


-2 


for M > 150 GeV. 


(B.4) 
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